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The physical properties of the Tomimatsu-Sato S — 2 spacetime are analyzed, with emphasis 
on the issues of the negative mass distribution in this spacetime and the origin of a massless ring 
singularity which are treated with the aid of an equatorially asymmetric two-body configuration 
arising within the framework of the analytically extended double-Kerr solution. As a by-product 
of this analysis it is proved analytically that the Kerr spacetime with negative mass always has a 
massless naked ring singularity off the symmetry axis accompanied by a region with closed timelike 
curves, and it is also pointed out that the Boyer-Lindquist coordinates in that case should be 
introduced in a different manner than in the case of the Kerr solution with positive mass. 
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The famous 6 = 2 Tomimatsu-Sato (TS2) solution [1, 2] discovered in 1972 within the framework of Ernst's complex 
potential formulation of the stationary axially symmetric gravitational field problem [3] describes the exterior field of 

■ a spinning mass different from that of a Kerr black hole [4] , and it has been discussed extensively in the literature. 
In the paper of Gibbons and Russell-Clark [5] it was suggested that the massive source of TS2 solution could be 

■ a specific rotating star; at the same time, the relevance of the solution to black-hole physics was discarded in that 
^ ■ paper because of the presence in TS2 spacetime of a region with closed timelike curves (CTC) containing a naked 

ring singularity outside the symmetry axis, and also because the hypersurface x = 1 was found not to be an event 
^ horizon of the solution. 

Several years ago, however, the interest in the global properties of TS2 spacetime in the black-hole context was 
revived by Kodama and Hikida [6] who pointed out that the poles a; = l,?; = ±lof the TS2 solution were in fact 
I two degenerate Killing horizons, whence they concluded in particular that the TS2 solution could provide "a new 
^ I possibility for the final states of gravitational collapse", the latter idea being recently reiterated in [7, 8]. In this 
' respect it is probably worth reminding that essentially all non-trivial members of the double-Kerr (DK) solution of 
Kramer and Neugebauer [9], to which the TS2 spacetime does belong, independently of their "pathologies" possess 
two Killing horizons (see, e.g., Ref. 10), p. 348), and, besides, that the binary black-hole configurations constitute 
[ only a 5-parameter subfamily of the general 8-parameter DK solution, to which the TS2 spacetime does not belong. 
Consequently, in our opinion, the most important physical aspects of the TS2 solution that are likely to be clarified 
first of all are the nature of the massless ring singularity and origin of the CTC region because providing a reliable 
explanation for these phenomena would also shed additional light on many other exact solutions with similar prop- 
erties. In the present paper we shall give a convincing demonstration that there exists a direct relation between the 
aforementioned two phenomena and the negative mass of the TS2 solution. This will be done, on the one hand, by 
obtaining analytical formulas for the mass distribution in the TS2 spacetime and, on the other hand, by considering 
^ . , the evolution of a one-parameter generalization of the TS2 solution, representing a specific non-linear superposition 
' of two Kerr-NUT sources, which will enable us to trace the emergence of the ring singularity on the stationary limit 
surface (SLS) of the constituent with negative mass. 

The paper is organized as follows. In Sec. 2 we give a brief review of the TS2 solution and some of its known 
properties. Here we observe in particular that the mass quadrupole moment of TS spacetimes has never been given 
correctly in the literature. In Sec. 3 the mass and angular momentum distributions of the TS2 solution are considered, 
establishing in this manner the presence of a negative Komar mass introduced by the non-zero NUT parameters of 
the corresponding DK solution. A one-parameter generalization of the TS2 solution which leaves the total mass and 
total angular momentum of the latter unchanged is constructed and analyzed in Sec. 4. This new binary configuration 
permits us, by varying the parameter zq from the larger to smaller values, to show that the appearance of the ring 
singularity on the SLS of the lower Kerr-NUT constituent takes place when the mass of this constituent changes its 
sign from plus to minus. The TS2 solution with negative total mass is considered in Sec. 5, and in this case two regions 
with CTCs appear, one touching the SLS outside the symmetry axis and the other not. Then the ring singularity 
arising only in the former region is subsequently given a generic interpretation as a locus of points at which the SLS 
touches the boundary of a CTC region. To support this definition of a ring singularity in the vacuum spacetimes we 
consider the Kerr solution as an additional example. Finally, Sec. 6 is devoted to discussion of the results obtained 
and conclusions. 
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II. THE TS2 METRIC AND ITS KNOWN PROPERTIES 

The TS2 metric in the Ernst-Perjes representation [11, 12] is defined by the Une element 



fidt-Ujdiff, (1) 



where the functions /, 7 and w are given by the formulas 



f = 



fx^ + fiu-il- 2/2) [(a;2 - 1)ct2 - ar] 



^ M^-(x2-l)(l-y2)^ 

^ n{l ~ y'')[(x'' - l)au + ^iT] 

- {x^ - - y^)a^ ' ^ 

and the polynomials /x, cr, v, r have the form 

= p-'ix^ - 1)2 + q^l - yY, <y = 2pq{x^ - y^), 
V = Ax{px^ +2x+p), T = -4gp"^(l - y2)(pa; + 1)^ (3) 

the real constants p and q being subjected to the constraint p^ -\- = 1. Prolate spheroidal coordinates x, y are 
related to the Weyl-Papapetrou cylindrical coordinates p, z by the equations 

I , . 1 



x=—{r++r-), y=—{r+-r-), r± = p"^ + {z ± n)'^ , (4) 

where k is the positive real constant in the line element (1). The inverse transformation is 

p = K\J{x'^ - 1)(1 - y2), z = Kxy. (5) 

The total mass M and total angular momentum J of TS2 solution are given by the expressions [1] 

2k- 

M = — , J = M^q, (6) 
P 

and, whereas the whole mass is contained in the central rod x — 1, y^ < 1, the angular momentum is distributed 
between the rod and the two poles a; = l,y = ±lin the following way [8]: 

7- 7 + 2 7 7 _ 2/^^(2 j _ i^'i^ + P) 

jPq ^ pq 

The poles were shown in [6] to be degenerate Killing horizons, which is in line with the general study of the DK 
solution made by Dietz and Hoenselaers [10]. The area of each horizon is [6] 

^pole = 87rKV^(l+p)- (8) 

The mass-quadrupole moment Q of TS2 solution deserves a special mention as it has never been published with a 
correct sign for nearly four decades. The easiest way to calculate Q is by using the expression of the corresponding 
Ernst potential, namely, 

£ = {A-B)/{A + B), 

A = p^x^ - 1) + g2(y4 _ 1) _ 2ipqxyix^ - y^), 

B = 2px{x^ -l) + 2iqy{y'^ -1). (9) 

Then the procedure [13] for the calculation of the Geroch-Hansen relativistic multipole moments [14, 15] readily yields 
the desired formula for Q: 

Q = -^M3(l + 3g2), (10) 



3 



and this Q turns out to be less than the respective quadrupole moment of the Kerr solution. Consequently, the 
rectified Q for all TS-solutions is 

Q = -M^ (^P^ + 9^) • (11) 

The geodesies in the TS2 spacetime were studied by various authors [2, 5, 7, 11], and it was precisely the investigation 
of the equatorial geodesies that turned out crucial for demonstrating that the hypersurface x = 1 was not an event 
horizon [5]. In the paper [16] it was shown that the singularities of TS2 solution defined by the invariants of the Weyl 
conform tensor coincide with zeros of the denominator ^ + B of the potential £ from (9). Therefore, here we can 
omit the explicit form of these invariants calculated in [16]. The function A + B takes zero values at the poles, x = 1, 
y = ±1, and at one point of the equatorial (y = 0) plane defined by the equation 

p^x^ + 2px{x^ - 1) - 1 = 0, (12) 

whence for the real root x = a;o > 1 we get the expression (p > 0) 




xo = ^[5-l + J3-d^--{l-2p^)], ^ = Jl + (2p)2/3(p2 - 1)1/3. (13) 



The point (xo, 0) determines the massless ring singularity outside the symmetry axis. The ergoregion of TS2 solution 
is bounded by two SLS defined by the equation / = 0, and the ring singularity lies on the inner surface if p > 0, i.e., 
if the total mass is a positive quantity. The existence of the region with CTCs was first pointed out in [5] , and this 
important region of causality violation is defined by the norm rj^ria of the axial Killing vector for which in [6] a rather 
cumbersome expression was given. However, a simple formula for rj'^rja was obtained long ago by Ernst [11], and we 
write it down in the explicit form, 

= ,^x^-i)ii-y^)r^-fu.^ = ^'^^-y'^^ , 

P^JJ 

N = p2(a;2 _ l)[p2(2.2 _ ^^^2 ^ 4^(p^2 +2x+p)+ q^{l - 2/2)2]2 

-V(l _ y^)[p^{x^ - l)(x2 - y^) + 2{px + 1)(1 - y2)]2^ (14) 

which is a good deal simpler than formula (34(i) of [6]. Note that the function D entering (14) is the denominator of 
the metric coefficient / already defined in (2). 

In Fig. 1 we have plotted a typical shape of the ergoregion, region with CTCs and location of the ring singularity (the 
point at which the curves determining the outer boundary of the CTC region and inner boundary of the ergoregion 
touch each other) for a particular choice of the parameters K,p > 0. 

Our idea to attribute the origin of the naked ring singularity and region with CTCs to the presence of a negative 
mass in the TS2 solution springs from the known analyses of the DK solution which provide a convincing evidence 
in favor of a direct connection between the negative mass of a rotating source and the corresponding ring singularity 
located on the SLS of that source [17-19]. Mention also that the recent paper [8] supplies one with an additional 
indirect hint on the negative mass in TS2 solution since the only reasonable interpretation of the fact that the poles 
x = I, y = ±1 have zero mass and non-zero angular momentum would be supposing that there are two sources 
associated with each pole which are endowed with equal positive and negative masses but unequal angular momenta. 
We are now turning to identifying analytically the negative mass distribution in the TS2 solution. 

III. THE NEGATIVE MASS DISTRIBUTION IN TS2 SOLUTION 

Since the metric function lu of the TS2 solution takes constant value on the interval {—k, k) of the symmetry axis, 
the distributions of the Komar mass and Komar angular momentum [20] on this rod can be analyzed with the aid of 
Tomimatsu's formulas [21] (both u) and to be evaluated a,t x = 1, y = z/k): 

M{z) = -^io[n{z)-n{-z)], J{z) = ^u{-2K-^w[n{z)-fl{-z)]}, (15) 

where z stands for the interval {—z,z), \z\ < k, of the symmetry axis containing the origin 2; = 0, and O is the 
imaginary part of the complex potential S of the TS2 solution: 

= -^, C=p^{x^-in2x''-y^ + l) + q'iy'-l)iy''-l). (16) 
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FIG. 1: The ergoregion (between two outer curves), the region with CTCs (between the axis and the inner curve), and the ring 
singularity in the TS2 solution (the particular choice of the parameters is k, — 2, p — 3/5, q — 4/5). 



Like in formulas (14), D is the denominator of the function /. 

On the interval |z| < k of the symmetry axis, uj and 51 take the form 



AKqz{z'^ 



pq (7^(z^ + K^)^ + 4k^(1 + p)^z^ ' 

so that for M(z) and J{z) we readily obtain 

Bk2(1 +p)z{z^ + K^) 



M{z) 



J3[g2(z2 ^ ^2)2 4^2(1 +p)2^2] ' 



j(,)^Mi±P)[M(z)-«] 



pq 



(17) 



(18) 



Restricting ourselves to the case of the positive total mass M, we must choose p > in the above formulas, and 
without lack of generality we can also put g > there as the negative values of q only change sign in the expression 
of the total angular momentum J. Then for the point z — z^ax at which M{z) takes its maximum value we easily 
find from (18) 



(19) 



SO that the corresponding value ^(zmax) yields the amount of the positive mass Afp too, while the amount of the 
negative mass Ad^ can be found either directly from (18) and (19) or simply by taking the difference between M 
and Afp (of course, this is justified because M{z) is monotonically increasing on the interval < z < Zmaxi a-nd 
monotonically decreasing on the interval z^ax < z < n). The resulting expressions for Afp and Af^ are 



A/p = — , Af„ = ^ ^, 

pq pq 



(20) 



while the angular momentum Jj-od of the rod is distributed between the momenta of the positive and negative masses 
in the following way: 



Jo 



2k2(2 - to) 

p2(l__p) 



Jn = - 



4k2(1 _ q) 
p^{l~p) ' 



(21) 
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Therefore, we have estabhshed analytically the presence of the negative mass in TS2 solution and identified its 
location on the symmetry axis (the segments p = 0, Zmax < \A < n)- The ratio \Mp/M^\ — 1/(1 — q) shows that the 
portion of the negative mass diminishes with increasing q, and tends to parity with the positive mass contribution in 
the limit g — > 0. 




FIG. 2: The mass distribution on the segment (— k, k) of the symmetry axis (the parameters are chosen as in Fig. 1). 




FIG. 3: The angular momentum distribution on the segment (— k, k) of the symmetry axis (the parameters are chosen as in 
Fig. 1). 



To illustrate the above general formulas, in Figs. 2 and 3 we have plotted the distributions M{z) and J{z) for a 
particular parameter choice k = 2, p = 3/5, q — 4/5 (the same as in Fig. 1), and one can see that on the interval 
(-2(2 - 73), 2(2 - V3)) the mass is positive, taking its maximum value Mp = 25/3, while the parts 2(2 — -\/3) < 
1^1 < 2 of the symmetry axis (the remaining parts of the rod) are characterized by a negative mass distribution with 
Mn = —5/3, the overall total mass of the rod being M — 20/3. 
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For the distribution of the angular momentum Jrod, shown in Fig. 3, we have Jp = —200/9, Jn = 760/9. To evaluate 
the total angular momentum, we also have to take into account the contribution of each pole Jpoie = ~40/3; then 

320 

J = Jp + Jn + 2 Jpole = -g- , (22) 

and it turns out that the contribution of the rotating negative mass into the angular momentum of the TS2 solution 
is the largest! 



IV. A ONE-PARAMETER GENERALIZATION OF THE TS2 SOLUTION 

Now our interest consists in finding a generalization of the TS2 solution such that it would contain an additional 

parameter by varying which wc could obtain evidence relating imambiguously the negative mass to the ring singularity. 
Mention that in the paper [22] Oohara and Sato showed how the TS2 solution can be obtained from the DK solution 
of Kramer and Neugebauer [9] in the special case of two equal corotating Kerr subextreme constituents. However, 
instead of considering the same equatorially symmetric binary configuration as in [22] we shall employ a different 
approach requiring the knowledge of the axis value of the potential (9) that will give us some additional useful 
physical information. 

We have found that on the upper part of the symmetry axis {y = 1, x = z/k) the potential £ of the TS2 solution 
(9) can be cast into the form {q > 0) 

This means, taking into account that the axis value of a single Kerr solution is given by the expression 

, ^ z — m — ia ,„ 
e{z) = — -, 24) 

z + m ~ la 

where m is the mass parameter and a the rotational parameter, that the TS2 solution can be formally considered 
as a non-linear superposition of two non-identical Kerr black holes - one with the mass ^(1 — ^/q) and angular 
momentum per unit mass —^^{1 — and the other with the mass ^(1 + ^/q) and angular momentum per unit 
mass ^v^(l + v^) - both black holes placed at the origin of coordinates and counterrotating. Since the terms (l±y^) 
are positive, then for p > the masses are positive too. Hence the question arises: where does the negative mass, the 
presence of which in the TS2 solution has been established earlier, come from? 

Probably an immediate simple answer to this question could be the following: a negative mass is needed to prevent 
two overlapping black holes from forming a unique Kerr black hole. On the other hand, it should be also remembered 
that we have a non-linear superposition in the presence of a strong gravitational field that can drastically change 
the initial interpretation of a binary configuration. To clarify the situation, we find it instructive to consider an 
immediate generalization of the axis data (23) of TS2 solution by introducing shifts of the black-hole constituents 
along the symmetry axis: 

^ ^-^o-g(l-Vg)(l-^Vg) + ^o-g(l + Vg)(l + i^) 
^^^> z-zo + ^{l-^){l + i^)' z + zo + !^{l + ^){l-i^y ^ ' 

where the first constituent is shifted from the origin at zq, and the second at — zq. Now we have a more general 
system of two non-identical separated Kerr black holes which in the limit zq leads to the TS2 configuration. 
A remarkable property of this system is that it has the same total mass and total angular momentum as the TS2 
solution, but its quadrupole moment generalizes formula (10): 

Q = -^M^il + 3q^) + Mzl (26) 

If we succeed to show that the solution defined by (25) contains an intermediate region with negative mass, then we 
may consider the TS2 spacetime as a final state of such special two-body system, with the negative mass redistributed 
within the rod and the poles in the limit zq = 0. So what we have to do now is to construct an exact solution 
corresponding to the axis data (25) and vary in it the parameter zq starting from some sufficiently large value 
and then go decreasing zq and calculating the Komar masses and angular momenta of the constituents and of the 
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intermediate region which separates them. As wiU be seen later, the lower constituent during such evolution will 
always remain subextreme, while the upper constituent is first subextreme, then becomes a hyperextreme object and 
eventually ends up as an extreme source at zq — 0. Therefore, we have to make use of the analytically extended 
double-Kerr (AEDK) solution [18, 23] obtained with the aid of Sibgatullin's integral method [24, 25] in order to be 
able to describe all possible changes of the binary system under consideration (the original DK solution is only valid 
for the description of the subextreme Kerr constituents). We recall that the complex potential £ of the AEDK solution 
has the form 

£ = E+/E^, 

E± = -(ai - Q!2)(a3 - ai){XiX2rir2 + X^X^rsr^) + (ai - a3){a2 - a^) 

x{XiX3rir3 + X2X4r2ri) - (ai - a4){a2 - a3){XiX4rir4 + X2X3r2r3), 
±{a3 - a4)[(ai - a3){ai - a4)X2r2 - {a2 - a3)(a2 - a4)Xiri] 
±{ai - a2)[(ai - a3){a2 - a3)X4r4 - (ai - a4){a2 - a4)X3r3], 

Xi = t~l'f:'fy r.^V^ + iz-a.r, (27) 
[at - Pi)[ai - P2) 

where /3; are arbitrary complex parameters, and a„ can take arbitrary real values or occur in complex conjugate pairs. 
Since Pi are simple poles in the axis data, then for our particular case we find from (25) 

A =^0- ^(l-^/^)(l + ^^/9), 132 = -zo - '^{1 + ,/q){l - i,/q). (28) 

On the other hand, the form of a„ corresponding to (25) can be obtained by solving (for z) the quartic equation 

(1 + q)[{z^ -K^- zlf - Ak'zI] + Sk'zo^z = 0, (29) 

which arises from the condition 



e{z) + e{z) = (30) 

of Sibgatullin's methos. 

Although Eq. (29) can be solved analytically, the resulting expressions are not wieldy to work with. Because we 
are basically interested in the qualitative picture of the evolution of the system, in what follows we shall assign some 
particular values to the parameters p, q and k, and consider different values of zq for which we shall be finding the 
corresponding concrete a„ by solving numerically Eq. (29). Once the values of Pi and a„ are known, the quantities 
Xn in (27) are also known, so that the corresponding metric functions and Komar characteristics of the binary system 
can be worked out with the aid of formulas (A2)-(A4) of Appendix. 

In table 1 we have described four particular configurations for the relatively large values of zo (the parameters p, 
q and k are fixed at 3/5, 4/5 and 1, respectively), where Mi, M2 and Mo denote the Komar masses of the upper 
constituent, lower constituent and intermediate region, respectively; Ji, J2 and Jo '^^'c the corresponding angular 
momenta. One can see that at first the two constituents are subextreme sources (all a„ are real) with positive masses, 
then, at smaller separations, the upper constituent becomes a hyperextreme object (a2 = (ii). The intermediate 
region has a small negative mass Mq introduced by the NUT parameters of the sources. No ring singularity appears 
in those configurations because the sources have positive masses. 

In table 2 the configurations undergo changes at smaller separations due to the spin-spin interaction of the three 
regions. The positive and negative masses of the lower constituent and intermediate region grow in value and then 
interchange their signs: at zq = 0.8, M2 is already negative, while Mq is positive. The ring singularity appears 
precisely after M2 becomes negative and is located on the SLS of the lower constituent. 

In Fig. 4 we illustrate the appearance and evolution of the ring singularity on SLS of the constituent with negative 
mass for different small zq, which leaves no doubt that the origin of this singularity is the negative mass. 

It should be especially remarked that the appearance of the region with CTCs in the configurations described by 
the above generalization of TS2 solution coincides with the appearance of the ring singularity and is also linked to 
the SLS of the lower constituent possessing a negative mass. We have checked numerically that the ring singularity 
always lies on the boundary of the CTC region, the latter boundary and the SLS of the lower constituent touching 
each other at the point determining the ring singularity. However, we have been unable to get high-quality plots of 
the particular CTC regions because of the corresponding complicated analytical expression of the norm of the axial 
Killing vector which was causing problems of the numerical computer processing. 
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TABLE I: Evolution of the generalized TS2 solution with varying zo («o > t)- The parameters p, q and /t are assigned the 
values 0.6, 0.8 and 1, respectively. 



Zq 


8 


6 


4 


2 


Mo 


-0.0015 


—0.0002 


—0.014 


—0.251 


Jn 


-0.258 


-0.385 


-0.719 


-2.824 


Ml 


0.168 


0.157 


0.139 


0.106 




—0.039 


—0.042 


-0.050 


—0.112 




^ 1 fi7 


^ 1 7fi 

O. J. 1 u 




^ /17Q 


J2 


9.186 


9.316 


9.658 


11.825 


/3i 


7.824 - 0.157i 


5.824 - 0.157i 


3.824 - 0.157i 


1.824 - 0.157i 




-11.157 + 2.824i 


-9.157 + 2. 824i 


-7.157 + 2. 824i 


-5.157 + 2. 824i 


ai 


8.048 


6.001 - 0.028j 


4.002 - 0.099i 


2.02 - 0.251i 


a2 


7.952 


6.001 + 0.028Z 


4.002 + 0.099i 


2.02 + 0.251i 


as 


-6.59 


-4.591 


-2.596 


-0.639 


a4 


-9.411 


-7.41 


-5.408 


-3.4 



TABLE II: Evolution of the generalized TS2 solution with varying zo {zo < k). The values of p, q, k are the same as in table 1. 





0.9 


0.8 


0.5 


0.001 


Mo 


-30.18 


18.115 


4.693 


3.334 


Jo 


-193.434 


112.614 


26.818 


15.664 


Ml 


0.084 


0.078 


0.055 


0.00016 


Ji 


-0.966 


-1.194 


-1.974 


-3.331 


M2 


33.430 


-14.86 


-1.414 


-0.001 


J2 


203.288 


-102.53 


-15.955 


-3.444 


/3i 


0.724 - 0.157i 


0.624 - 0.157i 


0.324 - 0.157i 


-0.175 - 0.157i 




-4.057 + 2.824i 


-3.957 + 2.824i 


-3.657 + 2.824i 


-3.158 + 2.824i 


Ql 


1.143 - 0.529i 


1.105 - 0.549i 


1.035 - 0.548i 


1 - 0.032i 


a2 


1.143 + 0.529i 


1.105 + 0.549i 


1.035 + 0.548i 


1 + 0.032i 




-0.01 


-0.039 


-0.222 


-0.968 


(Xa 


-2.277 


-2.172 


-1.847 


-1.032 



V. THE TS2 SOLUTION WITH NEGATIVE TOTAL MASS 

Once we have established that the TS2 solution with M > (p > 0) describes distributions of positive and negative 
Komar masses in which the positive mass prevails over the negative one, it seems quite logic to briefly discuss the 
opposite case when the negative mass prevails over the positive one and hence the total mass M becomes a negative 
quantity. It is interesting in particular to find out whether in the new situation we are going to have the same picture 
for the region with CTCs and the same location of the ring singularity as in Fig. 1. Like the previous case of positive 
M, the case M <Q {p <Q) can be readily treated with the aid of formulas (2), (3), (9), (12) and (14). 

In Fig. 5 we have plotted the SLS, the region of CTCs and the ring singularity of the TS2 solution for a particular 
parameter choice p = —3/5, q = 4/5, k = 2 which provides us with the same absolute value of mass M and the same 
angular momentum J as in the example used for plotting Fig. 1. It turns out that although the SLS in both figures are 
identical, the figures differ considerably in the shapes of the CTC regions and positions of the ring singularities. The 
first important distinction is the apcarance in Fig. 5 of the second region with CTCs which has a toroidal topology, 
and the second distinction is that the ring singularity in the case M < appears on the outer boundary of the 
ergoregion, thus reminding the position of ring singularities in the equilibrium configurations involving constituents 
with negative masses in the AEDK solution [18]. Observing that the first CTC region (having topology of a sphere) in 
Fig. 5 does not touch the inner SLS outside the axis and does not contain any ring singularity, we arrive at a generic 
conclusion that the massless ring singularity can be regarded as a direct outcome of the SLS touching the boundary of 
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the CTC region, and all this becomes possible only due to presence of a negative mass. Analytically, the singularity 
is defined in the prolate spheroidal coordinates by the point (a;o,0), Xq having the form 



a;o = < 




2p 



3-(52 + |(l-2p2)j , __L<p<o 



3-52- -(l-2p2) 



-l<p< 



1 

71 



(31) 



where 6 is the same as in (13). 

It looks interesting to compare Fig. 5 with an analogous figure in the case of a single Kerr solution with negative 
mass. An extra interest to do this consists in a necessity to correct a long-standing wrong belief that the Boyer- 
Lindquist coordinates [26] commonly used for the presentation and analysis of the Kerr metric are equally applicable 
to the situations with positive and negative masses. Then let us write down the metric functions /, 7 and co of the 
Kerr solution in the subextreme case using the Ernst representation in prolate spheroidal coordinates: 



f 



p^x^ + q^y^ - 1 

(px + 1)2 + g2y2 



=27 



p2(a;2 — 



2/tg(l- j/2)(pa; + i) 

p(p2a;2 -|- qi2y2 _ ' 



(32) 



the corresponding total mass and total angular momentum having the form 

K 



M = 



P 



(33) 
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FIG. 5: The regions with CTCs and the ring singularity in the TS2 solution with negative total mass. 

Like in the case of the TS2 solution, the negative values of M are defined by p < 0, and the Kerr solution for p < 
has a singularity at the point 



1 
P 



2/0=0 



Po 



J_ 

M 



Zq = 0, 



(34) 



as can be easily seen by considering the denominator of the function / in (32). The singularity lies off the symmetry 
axis (since > 1 or po > 0) in the equatorial plane and, therefore, is a naked ring singularity. It is also a curvature 
singularity, as follows from the form of the corresponding scalar invariant Ii of the Weyl tensor. 



/i = - 



2k,'^{px + 1 — iqy)^ ' 



(35) 



which can be calculated with the aid of the Economou-Ernst procedure [16] (the second invariant I2 is equal to zero 
identically). 

The SLSs of the Kerr solution with positive and negative masses are identical, being defined by the equation 



p^x^ + q^y^ -1 = 0, 



(36) 



and since the point x ~ — 1/p, y = fulfils (36), the ring singularity is located on the boundary of the ergoregion. 

On the other hand, the ring singularity also belongs to the boundary of the region with CTCs. Indeed, the norm 
of the axial Killing vector of the Kerr solution has the form 



N 



p'^Kpx + 1)^ + q'^y'^] ' 

[{px+ir+q^r-pV{x'~m-y% 



(37) 



and by a direct check one verifies that x = —1/p, y — satisfies N = identically, while the CTC region itself is 
defined by the inequality < 0. Hence, the ring singularity is a locus where the SLS and boundary of the CTC 
region touch each other. In the equatorial plane {y = 0) the region with CTCs extends from x = x^ to x — xi where 



1 



q 



xi = — l-A+:f- 



1 + — 92 



1/3 



(38) 



p \ 3A, 

and it can be readily verified that xi > xq and the norm rj^rja is negative definite on the interval xq < x < xi, y = 0. 
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The ring singularity of the Kerr solution with negative mass is massless, and all the mass is contained in the central 
rod j; = 1 of the symmetry axis, as can be readily seen with the aid of formulas (15) taking into account that the 
function ft of the Kerr solution is given by the expression 




2 4 6 

P 



FIG. 6: The SLS (black curve), region with CTCs (inside the red curve) and ring singularity in the Kerr solution with negative 
mass. 

In Fig. 6 we have plotted the SLS, the region with CTCs and the massless ring singularity (the point po = 8/3, zo = 0) 
of the Kerr solution by choosing the same particular values of parameters p = —3/5, q — 4/5, k = 2 as in the case of 
Fig. 5. Obviously, both figures 5 and 6 are quite similar qualitatively. 

Now, what happens if instead of analyzing the Kerr solution endowed with negative mass directly in the coordinates 
{x,y), one first passes to the standard Boyer-Lindquist coordinates (r, i9) via the formulas 

r- M , ^ 

x= , y — cosv, (40) 

K 

or simply starts with a representation of the Kerr metric in these coordinates? It is easy to see that while the choice 
M > yields an extension of the Kerr metric, the case M < leads to a contraction of the negative-mass Kerr 
spacetime because r = cos •& — then corresponds to location of a massless ring singularity which is a point off the 
symmetry axis! In this respect it is worth noting that Boyer and Lindquist in their paper [26] explicitly state that 
their consideration is restricted to the case M > 0, and in no place they mention that their coordinates are also a 
good option for M < 0. In any case, a reasonable coordinate change for M < might be (40) with M substituted by 
|M|, or by any other positive quantity which would have some physics behind it. Our own opinion is that the use of 
the spheroidal coordinates alone is sufficient for demonstrating the non-physical character of the Kerr solution with 
negative mass, especially taking into account that the singularity in this case is located outside the symmetry axis. 

The above said about the subextreme Kerr solution with negative mass can be summarized in the form of the 
following theorem. 

Theorem. The subextreme Kerr solution with negative mass develops a massless naked ring singularity off the 
symmetry axis defined by the point x = 1/p, y = (p — \ J/M\, z — in cylindrical coordinates) in the equatorial 
plane. The singularity is a locus where the SLS touches the region with CTCs, and the Kerr metric in this case does 
not describe a black hole. 
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VI. DISCUSSION AND CONCLUSIONS 



The analysis of the TS2 solution carried out in the present paper suggests that the physical interpretation proposed 
for this solution long ago by Gibbons and Russcll-Clark as describing the exterior field of a specific star whose 
prolateness is larger than that of a Kerr black hole still remains the most reasonable one. Indeed, the presence of 
a massless ring singularity, region with CTCs and negative mass distribution seems to discard the TS2 metric as a 
relevant candidate able to describe the final stage of the gravitational collapse; on the other hand, it is not excluded 
that a physically meaningful interior solution might exist which could be matched with the TS2 spacctimc on an 
appropriate boundary outside the problematic regions. 

While the analytical formulas obtained by us give a clear picture of how the negative mass is distributed within 
the central rod x = 1, < 1 of the TS2 solution, it is still not clear at the moment how to estimate the amount of 
negative mass contained in the poles x = I, y ~ ±1 except stating that it is equal in absolute value to the amount 
of positive mass. At the same time, it is not difficult to answer a question about where the negative mass comes 
from. The DK metric is a non-linear superposition of two Kerr-NUT solutions [27, 28], and the negative mass is 
introduced by the semi-infinite NUT sources, as it follows directly from the study of the latter sources either in the 
case of the pure NUT solution [29] or in the case of the combined Kerr-NUT spacetime [30] . The two NUT sources in 
the DK solution can be eliminated exclusively through imposing and solving the condition of asymptotic flatness and 
the axis condition. However, in the case of the TS2 solution only the former condition is satisfied, which means that 
some specific non-zero NUT sources are still present in this solution, giving rise to the negative mass. The TS2 axis 
data (23) well illustrates a simple idea that a product of two single Kerr data does not necessarily lead to a binary 
black-hole configuration, but rather to a specific system of two Kerr-NUT sources. 

We hope that we have been able to produce enough evidence in favor of an intrinsic connection existing between 
the negative mass, region with CTCs and massless ring singularity. The consideration of the TS2 and Kerr metrics 
with negative mass was crucial for establishing that it is precisely the 'contact' of the ergosphere and boundary of the 
region with CTCs that determines location of the naked ring singularity in the vacuum case. The fact that a correct 
interpretation of the negative-mass Kerr spacctimc has become possible only after a good understanding of a massless 
ring singularity of the TS2 metric demonstrates once again that the latter metric is a remarkable tool for testing new 
ideas in the field of exact solutions. 

Finally, we would like to remark that our discussion of the negative-mass Kerr solution casts doubts on the credibility 
of the claim made in the paper [31] about the existence of a negative-mass regular black hole surrounded by a massive 
fluid ring. Since the claim was made on the basis of a purely numerical approach to the problem and hence cannot 
be independently checked analytically, we restrict ourselves to only pointing out possible sources of error: (i) the use 
of the coordinates that do not allow to identify the presence of a ring singularity, (m) failure to reproduce unphysical 
regions by the numerical methods of paper [31], and (iii) automatic attributing the properties of event horizons to 
null hypersurfaces which can be not true [32]. In any case, a black hole with negative Komar mass, i.e., an object 
exhibiting repulsive gravity, looks to contradict the very essence of the notion 'black holes'. 
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Appendix A: Metric functions of the analytically extended DK solution 



The AEDK solution is the N = 2 specialization of the extended multi-soliton vacuum metric considered in the 
paper [23]. Its functions /, 7 and u) entering the line element 

Cis2 = /-I[e2')'(dp2 ^ ^^2-) ^ ^2^^2j _ ^^^^ _ ^^^-,2^ ^^^^ 

are defined by the following expressions: 

E+E_ + E+E_ 27 _ E+E_ + E+E_ _ 2i{GE_ - GE_) 
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G = -crA + (z + (j)r + ^ \ij{ai-\- ajYiTj -'^^i{ai' + Uji + ak'Yi 

l<i<j<4 i=l 

{i',j',k' j^i; i'<j'<k') 

4 



A 








l<i<j"<4 






= (-l)'+^(a,-aj)(«,. ^ 


- aj')XiXj, ^ i' < j') 




= (-l)Xai' - - 


ak'){af - ak')X„ {i' ,j' ,k' ^ i; 


a 


:= /3i+/32, A:= ^ 


Xij. 



(A2) 



l<i<j<4 



The explicit form of E±, Xi and is given in (27). 

The Komar mass and angular momentum of the subextreme constituent can be calculated with the aid of Tomi- 
matsu's formulas [21] 

Ma = [ cjfl^zdifdz, Ja = ^ [ w ( -1 - ^uiQ,;/] difdz, (A3) 



while for the calculation of the Komar characteristics of the hyperextreme object the more general formulas must be 
used [18]: 

Ma = ^(^ \p{\nf)^p-ujn^,\p^pjz + [p(ln/),^+w^p]^=^„dp 

Ja = "[2a;-2pw(ln/),p + (p2/"'+w')^.]p=Pot^^ 

- / [2pu;(ln/),, + (p2/"'+^')f^,p].=.„dp 
Jo 

rPo , 
+ [2pa;(ln/),, + (pV"'+w2)^p].=.,rfp), (A4) 

where Cl denotes the imaginary part of the Ernst potential £. For other notations we refer to papers[18, 21]. 



[1] A. Tomimatsu and H. Sato, Phys. Rev. Lett. 29, 1344 (1972). 
[2] A. Tomimatsu and H. Sato, Prog. Thcor. Phys. 50, 95 (1973). 
[3] F. J. Ernst, Phys. Rev. 167, 1175 (1968). 
[4] R. P. Kerr, Phys. Rev. Lett. 11, 237 (1963). 

[5] G. W. Gibbons and R. A. Russell-Clark, Phys. Rev. Lett. 30, 398 (1973). 
[6] H. Kodama and W. Hikida, Class. Quantum Grav. 20, 5121 (2003). 
[7] C. Bambi and N. Yosliida, Class. Quantum Grav. 27, 205006 (2010). 

[8] J. Gcgcnbcrg, H. Liu, S. S. Scahra and B. K. Tippctt, Class. Quantum Grav. 28, 085004 (2011). 

[9] D. Kramer and G. Ncugcbaucr, Phys. Lett. A 75, 259 (1980). 
[10] W. Dietz and C. Hoenselaers, Ann. Phys. (NY) 165, 319 (1985). 
[11] F. J. Ernst, J. Math. Phys. 17, 1091 (1976). 
[12] Z. Perjcs, J. Math. Phys. 30, 2197 (1989). 

[13] D. Fodor, C. Hoenselaers and Z. Perjcs, J. Math. Phys. 30, 2252 (1989). 

[14] R. J. Geroch, J. Math. Phys. 11, 2580 (1970). 

[15] R. O. Hansen, J. Math. Phys. 15, 46 (1974). 

[16] J. E. Economou and F. J. Ernst, J. Math. Phys. 17, 52 (1976). 

[17] C. Hoenselaers, Prog. Theor. Phys. 72, 761 (1984). 

[18] V. S. Manko V S, E. Ruiz and J. D. Sanabrisr Gomez, Class. Quantum Grav. 17, 1881 (2000). 

[19] J. A. Rucda, V. S. Manko, E. Ruiz and J. D. Sanabria-Gomez, Class. Quantum Grav. 22, 4887 (2005). 

[20] A. Komar, Phys. Rev. 113, 934 (1959). 



14 



[21] A. Tomimatsu, Prog. Thcor. Phys. 70, 385 (1983). 

[22] K. Oohara and H. Sato, Prog. Thcor. Phys. 65, 1891 (1981). 

[23] V. S. Manko and E. Ruiz, Class. Quantum Grav. 15, 2007 (1998). 

[24] N. R. SibgatuUin, Oscillations and Waves in Strong Gravitational and Electromagnetic Fields (Springer- Verlag, Berlin, 
1991). 

[25] V. S. Manko and N. R. SibgatuUin, Class. Quantum Grav. 10, 1383 (1993). 

[26] R. H. Boyer and R. W. Lindquist, J. Math. Phys. 8, 265 (1967). 

[27] E. Newman, L. Tamburino and T. Unti, J. Math. Phys. 4, 915 (1963). 

[28] M. Demiariski and E. Newman, Bull. Acad. Polon. Sci. Ser. Math. Astron. Phys. 14, 653 (1966). 

[29] V. S. Manko and E. Ruiz, Class. Quantum Grav. 22, 3555 (2005). 

[30] V. S. Manko, J. Martm and E. Ruiz, Class. Quantum Grav. 23, 4473 (2006). 

[31] M. Ansorg and D. Petroff, Class. Quantum Grav. 23, LSI (2006). 

[32] G. Gibbons and H. Kodama, Prog. Theor. Phys. 121, 1361 (2009). 



